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pH ' This talk discusses the relative impact of running-coupling and other higher-order corrections on the small- a; gluon- 

I gluon splitting function. Comments are made on similarities with some aspects of the Balitsky-Kovchegov equation, 

which arise because of the presence of an effective infrared cutoff in both cases. It is emphasised that, at least in the 
splitting-function case, the asymptotic small- a; behaviour has little relevance to the phcnomonologically interesting 
preasymptotic region. This is illustrated with the aid of a convolution of the resummed splitting function with a toy 
gluon distribution. 
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0^ ' 1 Introduction 

o 



I Detailed introductions to the more theoretical aspects of small- a; physics have been given in the contributions 
' to these proceedings by Ciafaloni [1] and by Mueller [2]. The former concentrated on our understanding of 
the all-orders perturbative structure of the linear problem of small-x parton multiplication, while the latter 
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' discussed the new phenomena that occur when the gluon density becomes so high that the small-x growth 
saturates. 

Oh, In the linear regime there have been extensive studies of the higher-order corrections. These are essential, 
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insofar as the leading-logarithmic (LLx) BFKL equation [3] for small- a; growth is clearly inconsistent with data 
(see for example [4,5]). However, the pure next-to-leading logarithmic (NLLa;) contributions to the evolution, 
[6,7], are so large that the problem appears perturbatively unstable. ° Techniques have been developed over 
the past few years to help understand the origin of the poor perturbative convergence, in the hope that one 
may then use that understanding to help reorganise the perturbative series into a more stable hierarchy. As 
discussed in [1], methods based on the combination of coUinear and small-x resummations [9-19] seem to be 
particularly successful in this respect. 

The situation in the context of saturation studies is less developed. Firstly, there is no definitive under- 
standing of how to extend linear LLx BFKL evolution to the saturation regime. One of the most widely studied 
models is the Balitsky-Kovchegov (BK) equation [20,21], which can be understood as resumming pomeron fan 
diagrams [22] , and for which an additional mean-field approximation is nearly always made (formally valid only 
for a thick nucleus, and over a limited energy range). Other more sophisticated approaches to saturation are 
currently being investigated {e.g. [23-26]), a number of which aim to account for pomeron loops, first shown 
to be important in some early numerical calculations [27,28] within the dipole approach [29]. Secondly, even 
within the simplest, BK, approach to saturation, studies of higher-order corrections have been less extensive 
than for the linear BFKL equation. 

One purpose of this talk is to examine some general lessons that have been learnt about the effects of 
higher-order corrections in the case of linear evolution and to discuss how they might be relevant also in the 
BK saturation case. 

A second part of this talk will examine briefly the outlook for phenomenological applications of the higher- 
order linear BFKL framework. 



*Talk presented at the QCD at cosmic energies workshop, Ettore Majorana Centre, Erice, Italy, September 2004. 
"Though there are certain observables for which specific implementations of the pure NLLx corrections may have reduced insta- 
bilities [8]. 
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2 General aspects of higher-order BFKL corrections 

In order to discuss higher-order corrections in hnear and saturating BFKL, it is important to understand what 
precisely to compare. A critical feature of the BK equation in this context is that its non-linear term provides 
an effective transverse infrared cutoff on the evolution, usually known as the saturation scale Qs- This cutoff 
scale varies as a function of y = Inl/a;. A very elegant formalisation of these properties has recently been 
given in [30]. 

Infrared cutoffs have long been investigated in linear BFKL. They arise (a) when one imposes them ad-hoc 
to eliminate the non-perturbative regime, or (b) implicitly, in the study of the gluon-gluon splitting function 
Pgg{z), which through factorisation contains just ultraviolet evolution, while infrared evolution (that below 
the factorisation scale) is entirely in the gluon distribution function, g{x,Q^), as illustrated in figure 1 (see 
e.g. [12,15,31]). 
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Figure 1. Evolution paths between two transverse scales, Qi and Q2, and their separation into the parton distribution (paths 
that go below Qi, blue if viewed in colour) and the perturbative DGLAP evolution (remaining paths, red). 
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In all BFKL-type problems the question of higher-order corrections is doubly complicated, because in 
addition to the usual NLLa; corrections (relative order as compared to LLx), the iteration of the kernel means 
that it is not possible to identify a unique scale at which to evaluate the kernel. In problems without cutoffs 
(and with two probes at similar hard scales), it turns out that it is nevertheless a reasonable approximation 
to use the hard scale of the problem as the effective scale, at least over a moderately large range of Y, 
because fluctuations in scale due to BFKL diffusion are, to a first approximation, symmetric around the 
hard scale [32-34]. In contrast, in each of the three cutoff-contexts mentioned above, it was discovered, 
independently [11,35,36], that if one uses the cutoff scale (the only physically unambiguously identifiable 
scale) as the renormalisation scale, then there are large negative corrections to the BFKL power, of relative 
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order as , which come about because the cutoff introduces an asymmetry: the BFKL evolution can only take 
place at scales larger than the cutoff, where the coupling is reduced by its running. 

That the correction should go as as can be seen as follows. Recall that in the usual fixed-coupling 
saddle-point approximation the gluon Green function between transverse scales k and ko at rapidity Y goes 
as 

G{Y;kM) - „^ , (1) 

where uj = asx(l/2), = OisNc/ir, and x, the BFKL characteristic function, and x" , its second derivative, are 
assumed to be evaluated at 7 = 1/2, unless otherwise stated. By examining dy \nG{Y, k,k) = cu — 1/2Y, one 
sees that the effective BFKL power receives a correction, Sui, of order 1/Y; the corresponding width in Infc^ 
of the solution is of order y/asX"Y- Equivalently if, for 'external' reasons, the width in In fc^ of the solution 
is limited to be At, then the BFKL evolution power will be suppressed by an amount 6uj ~ (a^x")/^*^ 
(c/. [37,38] for more precise calculations). When the evolution takes place with a running coupling, the width 
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is naturally limited by the cutoff in the infrared and by the low value of the coupling in the ultraviolet. The 
actual width of the solution is such that the two sources of suppression, i. e. the finite width of the solution and 
the running of the coupling, are of similar importance: (Q;sx")/At^ ~ wfoa^xAf, where b — {12Nc — 2n/)/127r. 
This gives At'^ ~ x"/(5q;sx)j or equivalently Suj ~ a^^^ {bx)'^^^x"^^^ ^ ^ correction of relative order ai^^ . 
This simple argument actually reproduces the whole of the suppression's leading functional dependence on , 
X, x" and b. 

Supplementing the above result with the relevant extra numerical coefficients, and additionally the NLL 
corrections, the small- a; power growth, ujc, of the Pgg splitting function at scale becomes 

tJc=i41n2a,(g2). - 4.0^2/3 - 6.5a, + O (a^/^) ) (2) 

One sees that, numerically, the running coupling and NLLx contributions are both large, negative, and of the 
same order of magnitude. In order to make a phenomenological prediction it is necessary to take into account 
the running of the coupling at all orders and to supplement the NLLx corrections with the yet higher-order 
coUinear-enhanced terms (which we refer to as NLLb). The results for the power are shown in figure 2. One 
sees that despite their different parametric dependence on as, in practise if one takes individually either the 
running coupling or the NLLb contributions, they lead to almost identical suppressions. Interestingly though, 
when taking both running and NLLb contributions, there is only limited extra suppression compared to either 
one individually.* 
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Figure 2. The small-x power growth lOc for the Pgg spHtting function in various approximations: LLx and NLLb (NLLx with 
additional enhanced higher-order corrections), each with fixed and running coupUng [13]. 



In the context of studies of the evolution of the saturation scale with both higher-order corrections and 
running coupling [40], a similar phenomenon has been observed, though there, the sequence of results that 
was shown corresponded to fixed-coupling LLx, running coupling LLx and running coupling NLLa; — if one 
considers just this combination of results then it is tempting (as was done in [2,40]) to make the statement 
that running coupling effects are dominant and the NLLx corrections are rather small. However, recently, the 
fixed-coupling (approximate coUinearly improved) NLLa; results were presented [41] for the evolution of the 

''An alternative way of viewing the results, [39], is to consider not the absolute change in lOc, but rather the fractional change as 
one includes various higher-order contributions — one then notices (with, say, as = 0.2) that the inclusion of a first higher-order 
contribution leads to about a 50% reduction in ujc, while the second has almost as large an effect, being a further 35% reduction. 
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saturation scale, and together with [40], those resuhs suggest that the picture is actuaUy very similar to the 
splitting-function case: each of NLLx and running coupling contributions are individually large and negative, 
but combining them leads to only a small amount of further suppression. 

Actually, such a result is quite natural: while at the lowest orders, e.g. Eq.(2), different sources of higher- 
order effects combine linearly, at higher orders there are strong non-linear effects. In the case of running 
coupling and NLLx effects there are actually three physical mechanisms at play: (a) since the cutoff causes 
the solution of the BFKL equation to be dominated by higher scales, where as is smaller due to its running, 
NLLx effects are reduced; (b) NLLx effects themselves reduce the dependence of lOc on a^, (suppressing it 
more at large as than at small as), slowing the running of cUc with transverse scale, as if there were a reduced 
'effective' /3-function, and this leads to a smaller running coupling correction; (c) the NLLx corrections cause 
a very strong suppression of the diffusion coefhcient, x", which means that limiting the width of diffusion, as 
happens due to the running of the coupling, has a smaller effect on the asymptotic power. 

The discussion so far has concentrated just on the power-growth of splitting functions and saturation 
scales. In the case of the splitting function, with the aid of recent technical developments, it has become 
possible to study the whole ^-dependence of the splitting function, even at preasymptotic values of x [13, 18]. 
Again, one can examine what happens when switching on, separately, running coupling and NLLa; effects, as 
shown in figure 3. As was the case when studying just the asymptotic power, lUc, one sees that, individually, 
running coupling and NLLx (or rather NLLb) effects are of similar magnitude. What should be noted here 
though, is that when considering the size of the (phcnomenologically relevant) preasymptotic region of x 
without growth, there is a rather large additional effect from the combination of running-coupling and NLLb 
contributions — for example the point at which the resummed splitting function starts to become larger than 
the LO DGLAP splitting function is a; ~ lO"-*^ for fixed-coupling LLx, x ^ 1Q~^ for running-coupling LLx or 
fixed-coupling NLLb, and x ~ 10^^ for running-coupling NLLb. 
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Figure 3. The Pgg splitting function with fixed and running coupling for the LLx and NLLb cases [13], compared to LO DGLAP. 

Figure 3 conveys what is perhaps one of the main general lessons to be retained from studies of resummed 
splitting functions: asymptotic properties of small- a; resummation have little relevance at today's energies. 
This statement holds in two senses: the behaviour of the splitting function at moderately small values of x is 
definitely not power-like; and general properties that one may deduce from studies of the asymptotic region 
{e.g. that combining NLLb and running-coupling effects provides only a modest extra suppression relative 
to each one individually) do not hold in the preasymptotic region. In the case of the splitting functions, the 
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specificity of tlie preasymptotic region can be traced to the appearance of new liierarchies in the perturbative 
structure, finite towers of terms a'si'^s In^ discussed in [1,42], leading to the characteristic dip structure 

at Inl/x- 1/^+0(1). 

In the case of the BK equation, a full study of preasymptotic effects including higher orders has yet to 
be carried out. It would presumably require that one know the structure of the NLLx terms not only for the 
linear part of the evolution, but also for the non-linear term."^ Currently however, the higher-order corrections 
to the non-linear term are not known. In the meantime it would nevertheless be of interest to have even just 
a full study of the x and structure of the BK-equation in which only the linear term was supplemented by 
higher-order corrections. It is to be noted though that some general information on the impact of higher-order 
corrections on preasymptotics in the BK-equation can already be obtained from studies [43-45] which solve 
the BK-equation in x, space with additional terms that partially mimic the linear NLLx corrections. 

3 Phenomenological impact of resummed splitting functions 

We have seen. Fig. 3, that preasymptotic effects are large in the resummation of the gluon-gluon splitting 
function, so much so that the BFKL growth only sets in for z ^ 10^^. This suggests that resummation may 
have only a modest impact on DGLAP fits. To determine robustly whether or not this is the case would 
however require that one carry out a complete DGLAP fit, with not only the Pgg splitting function, but also 
the whole matrix of splitting functions and the coefficient functions, preferably in the MS scheme, so as to aid 
comparison with existing fixed-order DGLAP fits, e.g. [46-48]. This represents a major programme of work, 
some aspects of which are currently being investigated. 




10"^ 10'^ 10"^ 10'^ 10"^ 10"^ 10"^ 10"'' 10' 



X 

Figure 4. Convolution of CTEQ6M [48] gluon with various splitting functions, normalised to the gluon itself, {Pgg ® g)/g. 

Nevertheless, some degree of insight into the possible phenomenological impact can be obtained simply by 
taking a fixed gluon distribution (here CTEQ6M [48], which has the advantage of being smooth at small x) and 
examining the convolution Pgg (g) g{x, Q^), shown in Fig. 4 normalised to g{x, Q^). As well as the convolution 
with the resummed (NLLb) splitting function, the plot shows the convolution with the fixed-order splitting 
function up to NNLO [49]. The comparison is to be used only for illustrative purposes since the fixed-order 
splitting functions are in the Ms scheme (though actually, at small x, the scheme is usually important only 



'^Whereas the universality features demonstrated in [30] ensure that the asymptotic properties of the solutions [40, 41] are inde- 
pendent of the details of the non-linear term. 
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starting from N'^LO), while the resummed splitting function is in the Qo scheme [50]. Furthermore at large x 
the NLLb resummation has been matched only to the LO DGLAP splitting function. 

In Fig. 4, because the gluon distribution itself rises at small-x, a feature of the splitting function at some 
given X value manifests itself in the convolution at somewhat smaller x. Thus, though the NLLb splitting 
function drops below the LO splitting function for x 10~^ (c/. Fig. 3, though the value there is different), 
this crossover in the convolution takes place at a; 10~^. For the crossover in the opposite direction the effect 
is much stronger, the NLLb splitting function overtaking the LO splitting function at a; ^ 5 • 10~^, whereas 
in the convolution this occurs below 10~^. 

Looking at the comparison with higher orders, one notices that at small x, the resummed convolution 
coincides quite closely with the NNLO convolution — this is perhaps not unsurprising, since down to a; ^ lO^"' 
there is a good deal of similarity between the NNLO and resummed splitting functions [1,42]. Only for 
x < 10^^ does one start to sec a difference between the NNLO and NLLb convolutions and, over the remaining 
phenomenologically accessible region, the NLLb convolution is intermediate between the NLO and NNLO 
results. If one is courageous {i.e. one believes that the main characteristics will remain the same after scheme 
changes and inclusion of the full matrix of splitting functions and the coefficient functions), one may take this 
to suggest that current NNLO fits [47] should be adequate down to a; 10~* and that only beyond does the 
fixed-order truncation truly start to break down.'' 

Acknowledgements. The splitting- function results described here have been obtained in collaboration 
with M. Ciafaloni, D. Colferai and A. Stasto. I am grateful to the Ettore Majorana Centre for the kind 
invitation and financial support to attend this very enjoyable workshop. 
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